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ABSTRACT 

Determinations of the solar oxygen content relying on the neutral forbidden transition at 630 nm depend 
upon the nickel abundance, due to a Ni I blend. Here we rederive the solar nickel abundance, using the same 
ab initio 3D hydrodynamic model of the solar photosphere employed in the recent revision of the abundances 
of C, N, O and other elements. Using 17 weak, unblended lines of Ni I together with the most accurate atomic 
and observational data available we find logeNi = 6.17 ± 0.02 (statistical) ± 0.05 (systematic), a downwards 
shift of 0.06-0.08 dex relative to previous ID-based abundances. We investigate the implications of the new 
nickel abundance for studies of the solar oxygen abundance based on the [O I] 630 nm line in the quiet Sun. 
Furthermore, we demonstrate that the oxygen abundance implied by the recent sunspot spectropolarimetric 
study of Centeno & Socas-Navarro needs to be revised downwards from logeo = 8.86 ± 0.07 to 8.71 ± 0.10. 
This revision is based on the new nickel abundance, application of the best available g/-value for the 630 nm 
forbidden oxygen line, and a more transparent treatment of CO formation. Determinations of the solar oxygen 
content relying on forbidden lines now appear to converge around logeo = 8.7. 

Subject headings: line: formation — line: profiles — Sun: abundances — Sun: atmosphere — Sun: photo- 
sphere — techniques: polarimetric 



1. INTRODUCTION 

The reference solar oxygen abundance has 
been revised over the past decad e from logeo = 
8.93 ± 0.04 (lAnders & Grevessd [l989h via 8.83 ± 0.06 
([Grevesse & Sauvall Il998[ IGS98I) to 8.66 ± 0.05 
(Usplund et al.l I2005L IAGS05I) . This downward slide 



has been brought on by the tandem influences of three- 
dimensional photospheric models, treatment of departures 
from local thermodynamic equilibrium (LTE), identification 
of blends, improved atomic data and better o bservations 
(lAllende Prieto et alJ l200l1: lAsplund et alJ 120041) . The new 
abundances of oxygen and other elements have solved 
many outstanding problems, but ruined agreement be- 
twe en helioseismologic al theory and observation (see 
e.g. iBasu & Antial (2008). This has prompted a reanaly- 
sis of photospheri c models, resulting in support for high 
dAvres et al.1 120061: ICenteno & Socas-Navarrol 120081: lAvresi 
20081) low (IScott et al.l 120061: ISocas-Navarro & Nortonl 
20071 iKoesteike et al.ll2008l: iMelendez & Asplundll2008l) and 
intermediate dCaffau et al .1120081) solar oxygen abundances. 

Many of these analyses rely upon the forbidden oxy- 
gen line at 630.0304 nm, know n to contain a significant 
blend from Ni I at 63 0.0341 nm (lAllende Prieto et alJl200H: 
iJohansson et alJ 120031) . The strength of this blend, and 
therefore the eo indicated by [Ol]630nm, depend crit- 
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ically upon the solar nickel abundance (em)- This is 
no less true of the ingen i ous sp ectropolarimetric work of 
ICenteno & Socas-Navarrol d2008l) than of any other study 
based on [Ol]630nm. Here we accurately redetermine eNi, 
and discuss the impact of the new value upon abundances 
from [O I] 630 nm. We show that eNi is model-d ependent, 
contradicting claims by ICenteno & Socas-Navarrol that their 
technique allows a nearly model-independent analysis of eo- 

2. MODEL ATMOSPHERES AND OBSERVATIONAL DATA 

We used the same 3D LTE model atmos phere and line 
format i on co de as in ear lier papers (e . g. lAsplund et alJ 
I2000bl 12004 . described by lAsplund eTatl OOOOat) . We per- 
formed comparative calcul ations wit h three ID models: HM 
dHolweger & Mulled fl974l) . marcs dGustafsson et al.lll975b 
lAsplund et al.l 1 19971) and 1DAV (a contraction of the 3D 
model into one dimension by averaging over surfaces of 
equal optical depth). Each ID model included a microtur- 
bulent velocity £ t = lkms -1 . We averaged simulated inten- 
sity profiles over the temporal and spatial extent of the model 
atmosphere, and compared results with th e Fourier Trans- 
form S pectrogr aph (FTS) disk -center atlas o flBrault & Neckel 
(fl987l see also [Neckel 1999J- We removed the solar gravi- 
tational redshift of 633 ms" 1 , and convolved simulated pro- 
files with an instrumental sine function of width Act = | = 
0.857 km s" 1 ^ reflecting the FTS resolving power R = 350000 
dNeckellll999l) . We obtained abundances with the 3D model 
from profile-fitting via a x 2 -analysis, fitting local continua in- 
dependently with nearby clear sections of the spectrum. For 
ID models we used the equivalent widths of 3D profile fits. 

3. ATOMIC DATA AND LINE SELECTION 

Our adopted Ni I lines and atomic data are given in TableQ] 
The paucity of good lines and atomic data in the optical pre- 
cludes any meaningful analysis of Ni II in the Sun. The most 
accurate Nil oscillator strength s come from the laborator y 
FTS branching fractions (BFs) of IWickliffe & Lawlerl (fl997h . 
put on an absolute scale with the time-resolved laser-induce d 
fluorescence (TRLIF) lifetimes of Bergeson & Lawlerl (fl993h . 
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TABLE 1 
List of neutral nickel lines 



Atomic levels Isotope Wavelength Ex. Pot. logg/ gf log7 ral j cr a W \ log en; Weight 



Lower 




Upper 






(nm, air) 


(eV) 


(eff.) 


ret. 








(pm) 


(3D) 




3d 8 ( 3 F)4s4p( 3 P) 


5 G 4 


3d 9 ( 2 D)4d 


3 G 5 




474.01658 


3.480 


-1.730 


WL 


7.899 


844 


0.281 


1.60 


6.18 


1 




r 1 




"0 


58 Ni 
60 Ni 


481.19772 
481.19926 


3.658 


-1.592 
-2.006 


J03 


8.285 






2. 12 


6.20 




3d 8 ( 3 F)4s4p( 3 P) 


5 G 2 


3d 8 4s( 4 F)5s 


5 F 3 




481.45979 


3.597 


-1.620 


WL 


8.053 


743 


0.236 


1.58 


6.17 




3d 8 ( 3 F)4s4p( 3 P) 


5 G 3 


3d 8 4s( 4 F)5s 


3 F 4 




487.47929 


3.543 


-1.450 


WL 


8.039 






2.35 


6.15 




3d 9 ( 2 D)4p 


3 D 2 


3d 8 4s( 4 F)5s 


5 F 2 




488.67108 


3.706 


-1.780 


WL 


8.211 






0.90 


6.13 


; 


3d 8 ( 3 F)4s4p( 3 P) 


5 G 4 


3d 8 4s( 4 F)5s 


5 F, 




490.09708 


3.480 


-1.670 


WL 


8.062 


693 


0.238 


1.79 


6.17 




3d 8 ( 3 F)4s4p( 3 P) 




3d 9 ( 2 D)4d 


3 G 4 




497.61348 


3.606 


-1.250 


WL 


7.962 


843 


0.282 


2.86 


6.13 


2 


3d 8 ( 3 F)4s4p( 3 P) 


3 F 4 


3d 8 4s( 4 F)5s 


5 F 5 




515.79805 


3.606 


-1.510 


WL 


8.093 


691 


0.236 


1.86 


6.13 


3 


3d 8 ( 3 F)4s4p( 3 P) 


G 5 


3d 8 4s( 4 F)5s 


St- 

F 4 




550.40945 


3.834 


-1.700 


WL 


8.063 


713 


0.240 


0.97 


6.18 


1 


3d 9 ( 2 D)4p 


l F 3 


3d 9 ( 2 D)4d 


3 G 4 




551.00092 


3.847 


-0.900 


WL 


8.215 






3.76 


6.17 


2 


3d 9 ( 2 D)4p 


'F 3 


3d 8 4s( 4 F)5s 


3 F 4 




553.71054 


3.847 


-2.200 


WL 


8.280 


695 


0.216 


0.31 


6.16 


3 


3d 8 ( 3 F)4s4p( 3 P) 


3 G 3 


3d 9 ( 2 D)4d 


3 G 4 


60 Ni 
58 Ni 


574.92795 
574.93039 


3.941 


-2.526 
-2.112 


WL 


7.944 


832 


0.284 


0.44 


6.17 


2 


3d 8 ( 3 F)4s4p( 3 P) 


3 F 4 


3d 9 ( 2 D)4d 


3 G, 


60 Ni 
58 Ni 


617.67980 
617.68200 


4.088 


-0.816 
-0.402 


WL 


8.162 


826 


0.284 


6.54 


6.17 


2 


3d 8 ( 3 F)4s4p( 3 P) 


3 F 4 


3d 8 4s( 4 F)5s 


3 F 4 




620.46048 


4.088 


-1.100 


WL 


8.244 


719 


0.247 


2.11 


6.19 


3 


3d 8 ( 3 F)4s4p( 3 P) 


3 F 3 


3d 9 ( 2 D)4d 


3 G 4 


60 Ni 

58 Ni 


622.39710 
622.39914 


4.105 


-1.466 
-1.052 


WL 


8.322 


827 


0.283 


2.79 


6.17 


3 


3d 8 ( 3 F)4s4p( 3 P) 


1' 


3d 9 ( 2 D)4d 


3 G 4 


60 Ni 
58 Ni 


637.82328 
637.82580 


4.154 


-1.386 
-0.972 


WL 


8.317 


825 


0.283 


3.20 


6.20 


3 


3d 8 ( 3 F)4s4p( 3 P) 


3 D 3 


3d 8 4s( 4 F)5s 


3 F 4 




641.45884 


4.154 


-1.180 


WL 


8.369 


721 


0.249 


1.68 


6.20 


2 



NOTE. — Wavelen gths and excitation potential s are fromjLitzen et al. 1 1993). Radiative damping is from VALD I Kupka et al. 1999). Transition 
designations are from Wickliffe & Lawler 1 1997}) ex cept in the case of 481. 2 nm, for which the designation is from VALD. References for gf- 
values are WL: IWickliffe & L awler ( 1997]) and J03: [Johansson et al. 12003). For lines with isotopic components g/'-values are effective only, 
rescaled to reflect the (terrestrial) isotopic fractions of Rosman & Taylor ( 1998). Collisional damping parameters a and a are courtesy of Paul 
Barklem (private communication, 1999; now in VALD. lBarklem et al.120001) . Equivalent widths are from profile fits using the 3D model. 



A s mall number of high-qu ality ^/-values are also available 
from lJohansson et"ail d2003l) . based upon FTS BFs and a sin- 
gle TRLIF lifetime. The uncertainties of the individual oscil- 
lator strengths we employ range from 0.02-0.07 dex, but most 
are accurate to ±0.03 dex. 

Nickel has five stable isotopes: 58 Ni, 60 Ni, 61 Ni, 62 Ni and 
64 Ni, present in the appro ximate ratio 74:28: 1 :4: 1 in the Earth 
dRosman & Tavlorlll998l) . Isotopic splitti ng of optical l ines is 
small, and dominated by 58 Ni and 60 Ni. ILitzen et all (1 1993b 
have obtained accurate laboratory FTS wavelengths for the 
isotopic components of many Ni I lines, which we include in 
TableQ]where applicable. We model such lines with two com- 
ponents, distributing the total oscillator strength according to 
the terrestrial 58 Ni: 60 Ni ratio. We are confident that the avail- 
able data sufficiently describe the isotopic broadening of solar 
lines, as laboratory FTS recordings are far better resolved than 
lines in the Sun. Isotopic structure has not been included in 
other determinations of the solar nickel abundance. 

We gave weightings to lines from 1 to 3 according to 
the absence of blends in the solar spectrum and the clarity 
of the surrounding continuum. Along with 16 unblended 
weak lines, we include the somewhat stronger Nil 617.7 nm 
(W\ = 6.54pm) owing to its clean line profile and its accurate 
atomic data. This line should be less affected by the rigors 
of strong line formation than single-component lines of the 
same equivalent width, thanks to the desaturating effects of 
isotopic broadening. To be conservative, we give this line a 
weighting of 2 . Our list has only two lines in common with 
iBiemont et al.l (119801) . mainly due to the absence of accurate 



^/-values for the 10 other lines used by those authors. 

We took wavel engths and excitation potentials from 
ILitzen etafl ( 199 3). and radiative damping from VALD 
dKupka et al.l 1 19991) . For most lines, we used collisional 
broadening parameters calculated for individual lines by 
Paul Barklem (priva te communication, 1999; now in VALD, 
iBarklem et al.l200"o|). For the remain der we interpolated in the 
tables of lAnstee & O'Maral (1 19951) and IBarklem & O'Maral 
d 19971) where possible . For Nil 488.7 nm we used the tra- 
ditional lUnsoldl (1 19551) formula with a scaling factor of 2.0 . 
Transition designations are from IWickliffe & Lawlerl [l99% . 
except for Nil481.2nm, where the designation is from 
VALD. Apart from the slightly stronger 617.7 nm line, all our 
lines are quite insensitive to the adopted collisional broaden- 
ing. The broadening treatment of Nil 488.7 nm thus has no 
impact on our abundance determination, nor does the ambi- 

uity in the identificat ion of the upper level of Ni I 48 1 .2 nm 

Johansson et al.ll2003l) . 

4. NICKEL RESULTS 

The mean nickel abundances we found using different 
model atmospheres are given in Table [2] Examples of pro- 
file fits to Ni I lines with the 3D model are given in Fig. Q] 
exhibiting similarly impressive agreement with observatio n 
as seen with other species (e.g. lAsplund et al.ll200"0al |2004|) . 
None of the models show abundance trends with equivalent 
width (Fig. 0, excitation potential nor wavelength, and the 
scatter is universally low, boding well for the internal con- 
sistency of all models. Very little difference exists between 



The solar Ni and O abundances 
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TABLE 2 

Logarithmic solar nickel abundances (mean ± 1 standard deviation) 



3D 


1DAV 


HM 


MARCS 


Meteoritic 


log e Ni (Ni I lines) 6.17±0.02 


6. 17 ±0.02 


6.26 ±0.02 


6.16±0.02 


6. 19 ±0.03 



X = 620.461 [nm] 
Ex. Pot- = 4.088 [eV] 
log e — 6.192 




620.46 
Wavelength [nm] 




622.38 


Wavele 


622.40 622.42 
ngth [nm] 




f ^\ 




/t : 


: Nil 






- \ = 637. S26 [nm] 






; Ex. Pot. = 4.154 [eV] 






- log € - 6,204 






637,80 


637.82 


637.84 



Wavelength [nm] 

FIG. 1. — Example spatially- and temporally-averaged, disk-center syn- 
thesized Ni I line profiles (diamonds), compared with observed FTS profiles 
(solid lines). We removed the solar gravitational redshift from the FTS spec- 
trum, convolved synthesized profiles with instrumental sine functions and 
fitted them in abundance and line shift. Arrows indicate the windows used 
for profile fits; in the uppermost panel, we used the entire synthesized region. 

3D and 1DAV abundances, implying that the mean tempera- 
ture structure rather than atmospheric inhomogeneities is the 
main reason for the difference between the 3D and HM re- 
sults. The 3D e Ni is in excellent agreement with the meteoritic 
value dAGS05l) . whereas the HM value is not. 

We adopt the 3D Ni I result as the best estimate of the solar 
abundance: 

e Ni = 6.17 ±0.05. 

The total error (±0.05 dex) is the sum in quadra- 
ture of the line-to-line scatter (±0.02 dex) and po- 
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FIG. 2. — Solar nickel abundances indicated by Nil lines, as computed 
with different model atmospheres. No substantial trends with line strength 
are visible with any of the four models. The similarity between 3D and 1DAV 
abundances suggests that the mean temperature structure of the 3D model is 
predominantly responsible for the difference between HM and 3D results, 
rather than the presence of temperature inhomogeneities. 

tential systemat ics arisi ng from the model atmosphere 
(±0.05 dex) . IAGS05I gave e N i = 6.23 ± 0.04, from 
iReddv et al.l (I2003I) using an ATLAS9 model (Kurucz, 
|http : / /kurucz . harvard . edu/ grids . html ). Pre- 
vious review s (e.g . IGS98I) adopted 6.25 ± 0.09, by 
IBiemont et all (119801) using the HM model. Our value is 
0.06-0.08 dex lower than earlier ones, and 0.09 dex less than 
our own HM-based estimate. There is presently no evidence 
for no n-LTE effects on our chosen lines in the Sun dAsplundl 
12005b . but without a dedicated study we cannot rule them out. 
After adjusting for ^/-values and equivalent widths, we find 
abundances 0.0 6 and 0.07 dex higher with the HM model than 
IBiemont et al.l for the two lines in common. We have not been 
able to trace the exact cause of these disparities, but tenta- 
tively attribute them to differences in radiative transfer codes, 
continuum opacities and implementations of the HM model. 

5. IMPLICATIONS FOR THE SOLAR OXYGEN ABUNDANCE 

The revised solar nickel abundance presented here has a di- 
rect impact upon any derivation of the oxygen abundance us- 
ing t he [O il 630 nm line, as this line is blended with one from 
Ni I (lAllende Prieto et al.ll200lT) . 
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ICenteno & Socas-Navarrol d2008l) used the Stokes V pro- 
file of [Ol]630nm to find an atomic ratio eo,atomic/eNi = 
210 ±24 in a sunspot. The y adopted an outdated gf -value 
for the [Ol]630nm line (cf. IStorev & Zeipp en 2000), caus- 
ing an overestimation of the ratio by 15% (+0.06 dex). They 
assumed logeNi = 6.23 to find eo,atomk and converted this 
to a bulk eo by calculating that 51% of oxygen resides in 
molecules. This is a reasonable assumption; in sunspots 
the only significant oxygen-bearing molecule is CO, which 
(roughly) forms as many molecules as there are carbon atoms 
available, due to the low temperatures. This number thus 
mirrors the assumed C/O ratio at the start of the calculation. 
That ratio only depends weakly on the choice of 3D or HM 
model, as seen in t he shift from 0.49 ±0.11 to 0.54 ±0.10 
between IGS98I and IAGS051 A more straightforward way of 
estimating the contribution from CO would be to say that 
the maximum eco is given by the adopted carbon abundance: 

e q ~ eo.atomic + £c- 

ICenteno & Socas-Navarrol claimed a nearly model- 
independent analysis because neither their CO correction 
nor nickel-to-atomic-oxygen ratio relied on an atmospheric 
model, and they believed cm to be well-established. The first 
statement is approximately true in the current debate, and the 
second is true of photospheric (but not sunspot) models. Here 
we have shown that cm is a mo del-dependent quantity, how - 
ever. The determination of eo via lCenteno & Socas-Navarrof s 
method is thus manifestly model-dependent, so there is no 
longer any reason to prefer placing a prior on the C/O 
ratio than on ec directly. Using our new nickel abunda nce, 
correcting the [O I] 630 nm gf, adopting the ec of lAGS05l and 
fully propagating all errors, we find an oxygen abundance 
of loge = 8.71 ±0.10 instead of their 8.86 ±0.07. Had 
we ad opted the traditional sunspot model of Malt by et al.l 
(1 19861) instead of the one in ferred from spectrum inversion by 
ICenteno & Socas-Navarrol we would have found 8.67 ±0.10. 
Retaining ICenteno & Socas-Navarrof s prior on the C/O ratio 
(with the error thereupon given by IAGS05I) . one would 
obtain 8 . 74 ±0.10 wit h their sunspot model and 8.66 ± 0.10 
with the iMaltbv et al. | model. Clearly their me thod is not as 
model-insensitive as lCenteno & Socas-Navarrol argued. 

Our new Ni abundance also m odifies analyses of 
|"0 1] 6 3 nm in the qu i et sola r sp ectrum. lAllende Prieto et all 
(l2001l) . lAsplund et al.l(l2004 and lAyresl(l2008l) all allowed the 



Ni contribution to var y freely in their 3D profile-fitting of the 
6 30 nm f eature, whilst lCaffau et al.l(l2008l) fixed it with the e Ni 
of lGS98l With the Ni abundan ces from Tableland th e labora- 
tory gf -value of the Ni I blend (iJohansson et al.l2003l) . we can 
now accurately predict the Ni contribution to [O I] 630 nm. In- 
dependent of the adopted ID or 3D model atmosphere, it is 
0. 17 pm in disk-center intensity, and 0. 19 pm in flux. In terms 
of oxygen abundance, this implie s a decrease by 0.04 d ex to 
logeo ~ 8.65 for the analysis of lAsplund et all (120041) . fur- 
ther improving the excellent agreem ent between d ifferent in- 
dicators. The derived abundance of lAvresI d2008l) wo uld de- 
crease to about 8.77 while that of ICaffau et all (120081) would 
increase to approximately 8.72. Because we now know the 
strength of the Ni blend, it is surprising that these two stud- 
ies yield different results for the remaining contribution from 
oxyge n, as they bo th rely on the same 3D C0 5 BOLD model. 
Since ICaffau et all employed several 3D snapshots whereas 
lAvresI used only one, we tentatively consider the former more 
reliable. No Ni abundance has yet been estimated with the 
CO s BOLD model, but regardless of its value our conclusions 
about the strength of the Ni I 630 nm blend, and thus its im- 
pact on oxygen abundances found by different authors, would 
remain uncha nged. The diffe r ence o f ap proximately . 07 dex 
in the revised lAsplund et ail d2004l) and ICaffau et alj (120081) 
abundances from [O I] 630 nm probably reflects the different 
mean temperature stratifications of the two 3D models. 

Given our reappraisal of th e oxygen abu n dance s of 
Centeno & Socas-Navarrol (120081) . lAsplund etaf] d2004l) and 



Caffauetal.1 (|2008|). toget her with the recent study of 



Melendez & Asplund (2008) using the [O I] 557.7 nm line, it 



now seems that results from forbidden oxygen lines are begin- 
ning to converge around logeo = 8.7. Whilst this agreement 
might come as a relief to some, it only serves to sharpen the 
current discrepancy between spectroscopy and helioseismol- 

ogy- 



We would like to take this opportunity to commemorate 
the work and life of Sveneric Johansson, his contribution to 
atomic spectroscopy in general and to nickel and [O I] 630 nm 
in particular. PS thanks IAU Commission 46, the ANU and 
the Australian Research Council for financial support. 
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